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A relativistic Coulomb-like resummation factor in QCD is 
suggested, based on the solution of the quasipotential equa- 
tion. 
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In describing a charged particle-antiparticle system 
near threshold, it is well known from QED that the 
so-called Coulomb resummation factor plays an impor- 
tant role Q. This resummation, performed on the 
basis of the nonrelativistic Schrodinger equation with 
the Coulomb potential V(r) = —a/r, leads to the 
Sommerfeld-Sakharov factor 



Sn 



1 - exp(-X nr ) 



7T a 



(1) 



which is related to the wave function of the continuous 
spectrum at the origin, |^(0)| 2 . Here v m is the velocity 
of the particle. An expansion of Eq. ([!]) in a power series 
in the coupling constant a reproduces the threshold sin- 
gularities of the Feynman diagrams in the form (a/v) n . 
However, in the threshold region one cannot truncate the 
perturbative series and the 5-factor should be taken into 
account in its entirety. A description of quark- antiquark 
systems near threshold, which has now been intensively 
investigated also requires this Coulomb resummation. 
The 5-factor appears in the parametrization of the imag- 
inary part of the quark current correlator, the Drell ratio 
R(s), which can be approximated in terms of the Bcthc- 
Salpeter (BS) amplitude of two charged particles xbs(x) 
at x = m . The nonrelativistic replacement of this am- 
plitude by the wave function which obeys the Schrodinger 
equation with the Coulomb potential, leads to the ap- 
proximation (P with a — » 4a s /3, for QCD. 

In the relativistic theory, especially for systems com- 
posed of quarks lighter than the top, the nonrelativistic 
approximation needs to be modified. To use the 5-factor 
within such a relativistic regime one usually uses the sim- 
ple substitution v m —> v with v = y/l — Am 2 / s. How- 



ever, the corresponding relativistic generalization of the 
5-factor is obviously not unique, for there are numerous 
ways of expressing the nonrelativistic velocity in terms of 
the relativistic energy y/s. For a systematic relativistic 
analysis of quark- antiquark systems, it is essential from 
the very beginning to have a relativistic generalization of 
the 5-factor. 

In this letter we suggest a new form for this rela- 
tivistic factor in the case of QCD. Our starting point is 
the quasipotential (QP) approach proposed by Logunov 
and Tavkhelidze H, in the form suggested by Kady- 
shevsky 0. To find an explicit form for the relativistic 
5-factor we will transform the QP equation from momen- 
tum space into relativistic configuration space The 
local Coulomb potential defined in this representation 
has a QCD-like behavior in momentum space Q. 

The possibility of using the QP approach to define the 
relativistic 5-factor is based on the fact that the BS am- 
plitude, which parameterizes the physical quantity R(s), 
is taken at x = 0, therefore, in particular, at relative time 
t = 0. The QP wave function is defined as the BS ampli- 
tude at r = 0, and the i?-ratio can be expressed through 
the QP wave function V'qp(p) by using the relation 



XBs(a; = 0) = J d£l p V>qp(p), 



(2) 



where dVt p = (dp) / [(2tt) 3 E p ] is the relativistic three- 
dimensional volume element in the Lobachevsky space 
realized on the hyperboloid E 2 — p 2 = m 2 .[] 

The QP equation in momentum space has the form 

(2E - 2E P ) V>(p) = j dQ k V (p(-)k) ^(k) . (3) 

The proper Lorentz transformation, Ak, means a trans- 
lation in the Lobachevsky space 



A k p = p(+)k = p + k 



Vl + P 2 



pk 



HVTTF 



(4) 



The role of the plane waves corresponding to the trans- 
lations (0) is played by the following functions 
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x In the following we will consider the case of two scalar par- 
ticles with the same masses m and use the system of units 
c = h = m = 1. 
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£(p,r) = (E p - pn) 



-1— zr 



(5) 



where r = nr and n 2 = 1. These functions correspond 
to the principal series of unitary representations of the 
Lorentz group and in the nonrelativistic limit (p < 1, 
r > 1) £(P7 r ) ~ * exp(ip • r). The orthogonality and 
completeness relations for these functions are 

J dfi p e(P,r)r(P,r')=*(r-r / ), 
| (dr) £(p, r) C (k, r) = (2^) 3 «S(p(-)k) , (6) 

where the relativistic momentum-space <5-function is 
5(p(-)k) = y/l + p 2 5(p - k). The QP wave functions 
in the momentum and relativistic configuration represen- 
tations are related as follows: 



(7) 



For a spherically symmetric potential the ^-transform 
of Eq. (H) is the equation 

J dn p (dr') (2E - 2E p ) C(p, v)C (p, r') ^(r') 
= V(r)ip(r), (8) 

where the right hand side is local. Here the transform 
of the potential is given in terms of the same relativistic 
plane wave, 



^(p)= f(dr)e (p,r)V(r) 



7(p(-)k)= /(dr)r(p(-)k,r)y(r). 



(9) 



The left hand side of this equation can be rewritten in a 
non-integral form by using the operator of the free Hamil- 
tonian [0 



- sinh (i— \ - ex (i — 
r \ dr I 2r 2 \ dr 



Hq = cosh I i — 
dr 



(10) 



where Ag jV is the angular part of the Laplacian opera- 
tor. The relation 7do£(p,r) = _E p £(p.r) allows one to 
re-express the equation in terms of finite differences 



(11) 



2E-2H )iP(r) = V(r)tp(r) 



This equation for the Coulomb potential has been in- 
vestigated in Ref. |10|. The solutions contain arbitrary 
functions of r with period i, the so-called the i-periodic 
constants, which appear in the solutions due to the fi- 
nite difference nature of the Hamiltonian (|l0|) . For some 
problems, such as defining the bound state spectrum, 
this i-periodic constant is not important. However, for 



the purpose of extracting the S'-factor, we must develop 
a method which avoids this ambiguity. 

Consider the Coulomb potential defined in relativistic 
configuration space 



V(r) = 



(12) 



The ^-transformation of Eq. (|T^) gives in momentum 
space the following function 



V(A) 



1 



Xa sinh xa 



(13) 



where the relative rapidity xa corresponds to A = 
p(— )k and is defined in terms of the square of the mo- 
mentum transfer by Q 2 = — (p — k) 2 = 2(cosIixa — !)• 
For large Q 2 the potential V(A) behaves as (Q 2 InQ 2 ) -1 , 
which reproduces the principal behavior of the QCD po- 
tential proportional to a s (Q 2 ) / Q 2 with a s (Q 2 ) being the 
QCD running coupling. 

According to Eqs. (§), (§), and (g), we find a rela- 
tion between the required BS amplitude and the QP 
wave function, xbs( x = 0) = ipQp(r = i). Performing 
a partial-wave analysis we further observe that the QP 
wave function for an ^-state will contain the generalized 
power (— r)^ +1 ) = i l+1 r(ir + l + l)/T(ir), which vanishes 
at r = i for £ ^ 0. Thus, we need only to consider the 
£ = wave function for which we can write ip( r ) = VK 7 ")- 
Introducing the function R(r) — rip(r) into Eq. (|^), we 
get 

- dx dr' sm X 'r smx'r'(2E -2cosh X ')R(r') 
t Jo Jo 

= V{r)R(r). (14) 

We will seek a solution of Eq. (|l4|) with the Coulomb 
potential ( |l2| ) in the form 



R(r) 



d(exp(ir()R((). 



(15) 



where the (^-integration is performed in the complex plane 
over a contour with endpoints (a, (3) |0]. Substituting 
Eq. (|T5|) into Eq. (Ill) we find the equationp] 



(16) 



\ d( exp(i r C) {2E-2 cosh Q R(Q 

J a 



dtexp(irQR(Q, 



which, when we integrate by parts, yields the two equa- 
tions 



2 We first perform the ^'-integration w ith a regularization 
factor exp(— ex' 2 ) an d then set e = after all calculations. 
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FIG. 1. Contour of integration in Eq. ( |T^ ) and singularities 
of the function (11 9) in the complex (,-plane. 



FIG. 2. Behavior of W-°'{s) with relativistic and nonrela- 
tivistic S-factors. 



exp(irC) (2E-2coshQ R(() 



C=/3 



(17) 



and 



(2E- 2 cosh C) R(0 =-aR((). (18) 



The solution of Eq. (18) is 



R(0 = C( X ) exp(C)[exp(C) - exp(- X ) 

1 -A-l 

exp(C) - exp(x) 



A-l 



(19) 



where A — ia/(2sinhx), E = coshx, and C(x) is an 
arbitrary function of \. 

The branch points of the function ( fl9| ) are ±x + 2nin 
(see Fig. |l|). The contour of integration must not inter- 
sect cuts which we take from — oo + 27rm to ±x + 2mn. 
In the case when the interaction vanishes, a — > 0, the so- 
lution R{r) should reproduce the known free wave func- 
tion sinx^/sinhx- Taking into account these remarks 
and Eq. ( |l7| ) for the boundary values at (a, /?), we take 
a = —R — is, (3 = —R + is with R — > oo. The vertical 
part of the contour to the right is given by Re £ = +R- 
It is also convenient for finding a connection to an in- 
tegral representation of the hypergeometric function to 
take the horizontal parts of the contour to be character- 
ized by Im£ = ±7r (see Fig. Q). 

The resulting solution does not contain the i-periodic 
constant and reads^] 



The representation of this solution in terms of the hyper- 
geometric function can be found, for instance, by the substi- 
tution x = x — In s. 



oo 

R(r) — C(x) sinh7rr / dx exp(irx) exp(x) 

-oo 
A-l 



X 



exp(x) + exp(-x) 



exp(x) + exp(x) 



-A-l 



(20) 



Comparing the asymptotic form of Eq. ( |20| ) at r — ► oo 
with the free wave function we can determine the con- 
stant C(x) and calculate |V>qp(«)| 2 which leads to the 
relativistic S- factor: 



x(x) 



1 - cxp l-X( X )} ' 



x(x) 



sinhx 



(21) 



where x is the rapidity which related to s by 2 cosh x = 
sfs. The function X(x) in Eq. (f2l"| ) can be expressed in 
terms of v as X(x) = ~kol\J\ — v 2 /v. 

We note that this new relativistic factor could have 
a significant impact in interpreting strong-interaction 
physics. In many physically interesting cases, R(s) oc- 
curs as a factor in an integrand, as, for example, for the 
case of inclusive r decay, for smearing quantities, and 
for the Adler D function. Here the behavior of S at in- 
termediate values of v becomes important. To illustrate 
the difference between the factors (|l|) and (21), in Fig. || 
we plot the principal contribution to R{s) for the vector 
currents, 



R«Hs) = V -^±S, 



(22) 



for the nonrelativistic case with w nr — > v and the rela- 
tivistic one, for a — 0.25. 

We conclude this letter by discussing the two limit- 
ing cases. In the nonrelativistic limit, v -C 1, the rela- 
tivistic S'-factor ( |2l] ) reproduces the nonrelativistic re- 
sult (|l|). In the ultrarelativistic limit, as it has been 
argued in Ref. |12], the bound state spectrum vanishes 
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as m — > because the particle mass is the only dimen- 
sional parameter. This feature reflects an essential differ- 
ence between potential models and quantum field theory, 
where an additional dimensional parameter appears. One 
can conclude that within a potential model, the S'-factor 
which corresponds to the continuous spectrum should go 
to unity in the limit m — > 0. Thus, the relativistic re- 
summation factor S obtained here reproduces both the 
expected nonrelativistic and ultrarelativistic limits and 
corresponds to a QCD-like Coulomb potential. 
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